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Abstract 

I We develop and validate a computational approach to nanomagnets. It is 

£SJ ' built on the spin wave approximation to a Heisenberg ferromagnet whose pa- 

rametcrs can be calculated from a first principles theory (e.g. density functional 
i theory). The method can be used for high throughput analysis of a variety of 

nanomagnetic materials. We compute the dependence of the magnetization of 
--.^ \ an iron nanomagnet on temperature, size and shape. The approach is applied 

to nanomagnets in the range of 432 atoms to 59 million atoms, a size which is 
^ ■ several orders of magnitude beyond the scalability of density functional theory. 

t3 ; 

1 Introduction 

O ■ 

Nanomagnets are promising materials for high density data storage devices, as they 
^ ' represent a bit by a single dot or magnetic domain. In such systems data storage 

densities greater than 10 12 bits/in 2 may be possible while being thermally stable at 
room temperature [1]. Important design choices for such systems are the material 
composition, the number of atoms, and the shape of the dot. All of these choices 
can strongly affect the magnetic properties of the system. However, an experimen- 
tal search of different combinations is time consuming. Thus, it is important to 
investigate efficient computational approaches to calculating magnetic properties. 
One might expect such computational approaches to reduce the cost of bringing 
nanomagnetic storage devices to market. 
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A significant effort is devoted to the calculation of magnetic properties of materi- 
als from first principles. These studies typically describe zero temperature properties 
of a small number of atoms (less than a thousand). A effective spin model with pa- 
rameters taken from such a first principle model can be simulated to predict the 
behavior of larger systems at higher temperatures. 

Thus from both fundamental and practical points of view, an efficient compu- 
tational approach to magnets is important. One such approach is given by the 
Heisenberg model with the Hamiltonian: 

H = ~^J(i - j)Si ■ Sj + gHBB cxt ■ J~]Sj + -ffcrystal-aniso, (1-1) 

i,3 i 

where Si is a three component spin at lattice site i = {ii,i2,h) satisfying the 0(3) 
invariant constraint 

Sf + Sf + Sf = S(S + l). 

In the above, g is the Lande g- factor, hb the Bohr magneton, B cxt an external 
magnetic field and S is the spin of the Heisenberg model. The first term is called 
the exchange Hamiltonian, the second the Zeeman Hamiltonian, and the third the 
crystal anisotropy. In this paper we shall restrict ourselves to the first term only. It 
provides a good description of the magnetization curve in zero magnetic field. We 
shall return to the physical effects of the remaining terms in a future study. The 
number of lattice sites or spins will be denoted by N = 2N x N y N z where N x , N y , N z 
are the number of atoms along a bcc lattice dimension for a rectangular paral- 
lelepiped. For the quantum Heisenberg model the spins do not commute and obey 
an algebra for which [S x , S y ] = iS z . 

The parameters J are exchange constants which can be calculated from a first 
principles model. Although the Heisenberg model can be formulated for any material 
where the J values are known, in this paper we study the quantum Heisenberg model 
with 5=1 used to describe bcc iron, which is an important material from a practical 
[2] as well as a fundamental point of view. The experimental data for the ratio of 
the magnetization of bcc iron to its value at T = is found in [3] to four digit 
accuracy. We will extend our analysis to other materials in future work. 

For bcc iron the J values have been calculated [4] and are listed in Table 1. Here 
one mRy corresponds to 13.6056923 meV and 157.887324 Kelvin. We plot these 
J values as a function of r = \i — j\ in Fig. 1. We can see that the couplings are 
ferromagnetic out to second nearest neighbor. 

Using somewhat different values for J, the Curie temperature was estimated 
to be 1A1AK in a mean field approximation and 950iT using the Random Phase 
approximation [5]. The latter is within 10 per cent of the observed value 1043-fT. 
An even better estimate, obtained using a local spin density approximation [6] gives 
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Table 1: J values from reference [4] 

1070K. 

The above results indicate that the Heisenberg model of ferromagnetism can be 
derived from first principles and gives a good description of ferromagnetism in bulk 
bcc iron. In the next two sections we apply this model to the description of bcc iron 
nanomagnets, which are of great experimental and theoretical interest. 

This paper is organized as follows. In Section 1 we summarize our computational 
method. In Section 2 we compute the finite volume effects in iron nanomagnets. In 
Section 3 we compute the finite shape effects in iron nanomagnets at fixed volume. 
In Section 4 we present our main conclusions. 

Our computational method can be summarized as follows: 

(1) Determine the parameters J of an effective Heisenberg model from a first 
principles zero temperature calculation and/or low temperature measurement. 

(2) Choose a lattice geometry describing the nanomagnet. This will determine 
its size and shape. In this paper we will choose bcc iron and a parallelepiped lattice 
geometry, but other choices are possible. 

(3) Apply the Holstein-Primakoff transformation to map the spin variables to a 
set of harmonic oscillators whose nearest neighbor interactions determine a lattice 
Laplacian. 

(4) Solve for the eigenstates and eigenvalues of this lattice Laplacian. 

(5) Use these eigenvalues to determine the magnon dispersion relation. 

(6) The main computational step is to construct a sum over these eigenvalues 
using Bose-Einstein statistical mechanics to determine the magnetization, given by 
the expectation value M = (J2 Sf) as a function of temperature. We define the 

i 

thermal average by (O) = tr (Oe~ H / kT )^j for an operator O. 

(7) Vary the temperature, size and shape of the nanomagnet to determine its 
range of magnetic properties. 
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Figure 1: J values from [4] plotted as a function of r = \i — j\. 

(8) Compare the above derived magnetization curves against experiment. 

We comment on the main steps (3) and (4). Since we consider only nearest 
neighbor terms, the Holstein-Primakoff transformation, which is a linearization of 
the full many body Hamiltonian, yields nearest neighbor coupling only, and a dis- 
crete version of the Laplace operator. The standard continuum Laplacian is an ap- 
proximation which ignores atomic spacing and lattice geometry. The construction 
thus is applicable to arbitrary lattices, including lattices with defects and arbitrary 
shapes. For such problems, the eigenvalues of the Hamiltonian are solved numer- 
ically. For simplicity we consider here a rectangular geometry and an absence of 
defects. In this case the eigenvalues are known in closed form. 

A similar approach can be found in the description of spin-waves in ribbon 
shaped iron nanoparticles in [7]. An earlier study of the finite size effects, based 
on the analysis of the spin-wave spectrum of Heisenberg spins was performed in 
[8]. The inclusion of higher order spin wave interactions can also be considered 
computationally [9]. These become more important as one approaches the Curie 
temperature [10] and can included in full quantum monte carlo simulations of the 
quantum Heisenberg model [11]. 
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2 Finite size effects in bcc iron nanomagnets 

The presence of the T 3 / 2 term in the empirical formulas for bcc iron indicate that 
quantized spin-waves or magnons and the associated Bloch law are physically very 
important [12, 13]. 

The magnon or spin-wave appoach is most easily introduced using the trans- 
formation of spin variables of Holstein and Primakoff [14] . This transformation is 
given by 

S+ = V2Sfidi 
S~ = V2Saffi 
Sf = S - afcii 

where aj and af satisfy the usual harmonic oscillator algebra and fi, Sf, and 
are defined by: 

fi = (1 - afa^S) 1 / 2 
S+ = Sf + iSf 

q— — qx a qV 

In terms of these variables and neglecting terms higher order than quadratic the 
Hamiltonian becomes 

H = ~ E J (* " 3) ~ 2 - i){afa 3 - a+o,) + 0(a 4 ) 

Treating the above form of the Hamiltonian as defining a lattice Laplacian one solves 
for the eigenfunctions and eigenvalues through the equation: 

J2((i\H\j) -^)Cj = 
j 

The simplest case is the free or discrete-Neumann boundary condition where the 
spin on the boundary can rotate freely. In this case the eigenfunctions are given by 

Cf = ^cos(^(n - i))cos($-(i 2 - i))cos(^(z 3 - i)) 



where 



and 



h = 0,1, 2,..., 2N X -1 
£ 2 = 0,1, 2,..., 2N y -l 
4 = 0, 1,2,. ..,2N Z -1 



m K 2N X (1 + S tl ,o) ' y 2N y {l + 5, 2fi Y K 2N z {l + 5 hfl )> 
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The eigenvalues are given by: 



o t fa a nil ^ 2 ^ 3 ^ foi\ 
uj£ = 2.J 8 — 8 cos — — cos — — cos — — (2.1) 

V 2N * 2N v 2Nz J 

The eigenfuction analysis is similar to [15] except that in that reference discrete 
Dirichlet boundary conditions and eigenfunctions given by sines were used. 

Given these basis functions one can transform the a« variables to a new variable 
be via: 

k 

In terms of the tranformed variables the Hamiltonian becomes: 

H = Y J ^th + o(b A ) + ..., 

and the magnetization is given by: 

- 1 - WS £ <W - 1 - WS £ exp( ,,/T) - 1 - 

The Bloch T2 law for bulk iron follows from treating i values as continuous with 
magnitude from zero to infinity. For sufficiently small magnets, the compact and 
discrete aspects of the magnon dispersion relation and finite size modifications to 
the Bloch law can be observed. 

Given the form of the magnetization it is useful to define an intrinic temperature 
for the bcc iron system given by T a = D/a 2 = 2J/kb = 395.902i"C where D is the 
spin stiffness coeficient given by D = .281eVA [7]. The lattice spacing for bcc iron 
is a = 2.87 A = .287nm. 

In Fig. 2 we plot the quantity M(T)/M(0) using (2.1) for different size magnets. 
The basic effect is that the magnetization increases for a given temperature if we 
decrease the size of the nanomagnet. The behavior is similar to the magnetization 
of nanoscale iron islands studied in [16]. The finite size effect is larger as one goes 
to higher temperatures. Beyond 400K the free magnon approximation breaks down 
and one has to include magnon interactions. 

For nanomagnets of size greater than (28.7nm) 3 or 100 atoms on a side the 
finite size effects are small and one can treat them using a 1/N X expansion. For 
the case N x = N y = N z the parameter 1/N X completely determines the size of the 
nanomagnet. It is clear from Fig. 3 that there is a decrease in magnetization per 
atom as the size increases. For large values of 1/N X the curve can be empirically fit 
to the formula: 
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Figure 2: Finite size effects in iron nanomagnets. Calculated M(T)/M(0) for bcc 
iron nanomagnets using (2.1). The nanomagnets are of size (1.722nm) 3 with 432 
atoms and (2.87nm) 3 with 2000 atoms. The experimental magnetization curve for 
bulk iron [3] is shown for reference. 



= 1 - .031245(T/T a ) 3 / 2 + .056431(T/T a )^ + O(^) (2.3) 

Note as M(0) is proportional to volume, the first and second terms in the above fit 
can be interpreted as volume and surface area contributions to the magnetization. 
The curve is nearly linear versus 1/N X in the range between N x = 600 and N x = 100. 
We also studied the behavior between N x = 100 and N x = 10 (see Fig. 4) and also 
found that the same linear dependence is valid. This disagrees with the analysis of 
[7] [4th formula in Appendix B] who found that the leading finite size correction 
was proportional to (T/T a ) 5//4 ]y 2 . A possible source of error in the analysis of this 
correction is due to the use of the integral transformation 

1 
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Figure 3: Calculated M(T = 295K)/M(0) for bee iron nanomagnets as a function 
of 1/N X , the reciprocal of the number of atoms on a side under the condition N x = 
N y = N z for the range N x = 600 to N x = 100. 
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Figure 4: Calculated M(T = 295K)/M(0) for bcc iron nanomagnets as a function 
of 1/N X , the reciprocal of the number of atoms on a side under the condition N x = 
N y = N z for the range N x = 100 to N x = 10. 
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For N large one should approach the magnetization curve of bulk iron. In Fig. 
5 we plot the magnetization calculated using the magnon approximation for an iron 
nanomagnet with 639 x 639 x 639 atoms on a side, with volume (183nm) 3 and 
521,834,238 atoms. It is clear from the figure that additional effects such as spin 
wave interactions are physically important beyond temperatures of 400K. 
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Figure 5: Magnon calculation of the magnetization for bcc iron nanomagnets with 
639 x 639 x 639 atoms on a side, with volume (183nm) 3 and approximately 522 
million atoms compared to bulk experimental curve of [3]. 



3 Shape effects in bcc iron nanomagnets 

One of the defining characteristics of a nanomaterial is that nanomaterial properties 
differ from those of the bulk. In the previous section we studied the finite size effect 
on the magnetization at fixed shape. In this section we fix the volume and study the 
effect of different shapes on the magnetization. Although we only consider shapes 
of the form of a parallelepiped with N x x N y x N z atoms on a side, it is possible to 
consider more complicated shapes by modifying the discrete Laplacian which leads 
to the magnon dispersion relation. 

For the simple case when two of the lengths of the parallelepiped are equal say 
N x = N y there is only one independent shape parameter which we take to be h = 
N z /N x using the notation of [7]. At fixed volume the magnetization only depends on 
this parameter and temperature. In Figure 6 we plot the magnetization as a function 
of temperature for various shapes under the assumption of N x = N y . The dramatic 
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dropoff in magnetization of the nanowire is in accord with the expectation of the 
absence of ferromagnetism in one and two dimensions [17] regulated by the finite size 
of the magnet. An infinite one dimensional magnet would have zero magnetization 
at any finite temperature. Indeed the one dimensional quantum Heisenberg model 
was exactly solved in [18] and exhibits an absence of ferromagnetism. 

To further illustrate this phenomenon we plot the magnetization in Fig. 7 
at 100K and 295K as a function of the shape parameter h with total volume 
(34.44nm) 3 . The data can be fit to the equation: 

= 1 - .031245(T/T a ) 3 / 2 + .056431(T/T a )^ 

+.03831 {T/T a ) 1 ^(2h 1 / 3 + h- 2 / 3 - 3) - .02083(T/T a )^(2/i- 2 / 3 + h A l 3 - 3) 

(3.1) 

Here we have defined N m = (N x N y N z ) 1 / 3 . This equation reduces to (2.3) for N x = 
N y = N z and h = 1. 

To gain further understanding of the various terms in this equation consider the 
more general case of a nanoribbon where two lengths need not be equal and two 
aspect ratios are defined by h\ = N z /N x and hi = N z /N y . In Figure 8 we plot 
nearly 2000 different shapes and fit to the equation: 



S = 1 - -031245(T/T a ) 3 / 2 + .056431(T/T a )£ 

% 1/3 + hT h- 1/3 + h- 1/3 h 

-.02083{T/T a ) 1 f(h~' i/3 h 2 2 /3 + h~ 4/3 hl /3 + h 2 / 3 h 2 2 /3 - 3) 



+.03831(T/T a ) 7 ^(/ i 2/3 /i 2 1/3 + h 2 2 /3 h' 1/3 + h~ 1/3 h~ 1/3 - 3) (3.2) 



Although complicated, this formula can be understood by considering an expansion 
of the magnetization before division by M(0). Defining L x = aN x , L y = aN y , L z = 
aN z and (5 = 1/kT, the second term on the right side is of the form L x L y L z / (Df3) 3 / 2 
and is responsible for the usual T 3 / 2 law. The remaining terms are linear combina- 
tions of three types: 1 = (L x L y L z ) 2 / 3 /(D0), 2 = {L x L y + L y L z + L z L x )/(D0) 
and O3 = (L 2 + L 2 + L 2 )/ (D{3) . These are the unique terms which are quadratic in 
the L's and invariant under permutations of the axes of the parallelepiped. Express- 
ing the L's in terms of the parameters N m , hi and hi and dividing by M(0) (which 
is proportional to N x N y N z ) we obtain an expression of the form (3.2). Equation 
(3.2) reduces to (3.1) when h\ = hi- As before (3.1) reduces to (2.1) when h = I. 

An experimental determination of the magnetization of iron nanoribbons was 
carried out in [7]. In that work the long direction of the nanoribbon was not specified. 
As in that work we use the spin-wave model to compare with this data and fit the 
long direction of the nanomagnet. The prediction of the spinwave model for the 
nanoriboon is shown in Fig. 9. Our determination of the length of the long direction 
is 3645 nm. This is of the same order of magnitude as the 2000 nm determination in 
[7] where the spin wave method, but not the bcc lattice finite size dispersion relation 
was used. 
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Figure 6: Shape effects in nanomagnets. M(T)/M(0) for bcc iron nanomagnets of 
different shape, each consisting of 3,456,000 atoms with total volume (34.44nm) 3 . 
The shape from top to bottom are 120 x 120 x 120; 20 x 20 x 4320; 15 x 15 x 7, 680; 
12 x 12 x 12, 000 and 10 x 10 x 17, 280 atoms on a side. 

4 Conclusion and Future Direction 

We developed and validated an efficient computational approach to nanomagnets. 
The method takes input from a first principle calculation of the Heisenberg model 
parameters and uses a magnon approximation to the quantum Heisenberg model to 
calculate the magnetic properties. Although we demonstrated the method for iron, 
this type of calculation can be carried out for other materials and remains valid at 
room temperature. Composite systems such as iron-rhodium can be treated by our 
method as first principle calculations of the J values are known [19] in that case. 
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